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We consider the interaction potential between two chiral rod-like colloids which consist of a thin 
cylindrical backbone decorated with a helical charge distribution on the cylinder surface. For suf- 
ficiently slender coiled rods a simple scaling expression is derived which links the chiral 'twisting' 
potential to the intrinsic properties of the particles such as the coil pitch, charge density and electro- 
static screening parameter. To predict the behavior of the macroscopic cholesteric pitch we invoke 
a simple second- virial theory generalized for weakly twisted director fields. While the handedness 
of the cholesteric phase for weakly coiled rods is always commensurate with that of the internal 
coil, more strongly coiled rods display cholesteric order with opposite handedness. The correlation 
between the symmetry of the microscopic helix and the macroscopic cholesteric director field is 
quantified in detail. Mixing helices with sufficiently disparate lengths and coil pitches gives rise to a 
demixing of the uniform cholesteric phase into two fractions with a different macroscopic pitch. Our 
findings are consistent with experimental results and could be helpful in interpreting experimental 
observations in systems of cellulose and chitin microfibers, DNA and fd virus rods. 

PACS numbers: 61.30.Cz, 64.70.M, 82.70.Dd 



I. INTRODUCTION 

In contrast to a common ncmatic phase, where the ne- 
matic director is homogeneous throughout the system, 
the cholesteric (chiral nematic) phase is characterized by 
a helical arrangement of the director field along a com- 
mon pitch axis. As a result, the cholesteric phase pos- 
sesses an additional mesoscopic length scale, commonly 
referred to as the 'cholesteric pitch which characterizes 
the distance along the pitch axis over which the local di- 
rector makes a full turn [l[ . The behavior of the pitch as a 
function of density, temperature and solvent conditions 
is of great fundamental and practical importance since 
the unique rheological, electrical and optical properties 
of cholesteric materials are determined in large part by 
the topology of the nematic director field. 

In recent years considerable research effort has been 
devoted to studying chirality in lyotropic liquid crystals 
which consist of colloidal particles or stiff polymers im- 
mersed in a solvent. In addition to a number of synthetic 
helicalpolymers such as polyisocyanates [l,!!] and polysi- 
lanes (j] which form cholesteric phases in organic solvents 
there is a large class of helical bio-polymers which are 
known to form cholesteric phases in water. Examples are 
DNA 0, @ and the rod-like fd -virus 0, polypeptides 
d, II] , chiral micelles [l^j , polysaccharide s flU , and mi- 
crofibrillar cellulose (and its derivatives) [12| and chitin 
[l3| . In these systems, the cholesteric pitch is strongly 
dependent upon the particle concentration, temperature 
as well as e.g. the ionic strength which has been the 
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subject of intense experimental research [sl. [l3 - l23j . 

Understanding the connection between the molecular 
interactions responsible for chirality on the microscopic 
scale and the structure of the macroscopic cholesteric 
phase has been a long-standing challenge in the physics 
of liquid crystals [H . The chiral nature of most biomacro- 
molecules originates from a spatially non-uniform dis- 
tribution of charges and dipole moments residing on 
the molecule. The most prominent example is the 
double-helix backbone structure of the phosphate groups 
in DNA. Combining the electrostatic interactions with 
the intrinsic conformation of the molecule allows for a 
coarse-grained description in terms of an effective chi- 
ral shape. Examples are bent-core or banana-shaped 
molecules (23 . [25| where the mesogen shape is primarily 
responsible for chirality. Many other helical bio-polymers 
and microfibrillar assemblies of chiral molecules (such as 
cellulose) can be mapped onto effective chiral objects 
such as a threaded cylinder |2^ . [26l [27| . twisted rod 
[H, Hi] or semi-fiexible helix 

The construction of a microscopic theory for the 
cholesteric phase is a serious challenge owing to the com- 
plexity of the underlying chiral interaction and the in- 
homogeneous and anisotropic nature of the phase [30j . 
Course-grained model potentials aimed at capturing the 
essentials of the complex molecular nature of the electro- 
statics of the surface of such macromolecules have been 
devised mainly for DNA [sH - fssj . A more general elec- 
trostatic model potential for chiral interactions was pro- 
posed much earlier by Goossens [s^ based on a spatial 
arrangement of dipole-dipole and dipole-quadrupole in- 
teractions which can be cast into a multipole expansion 
in terms of tractable pseudo-scalar potentials j37| . This 
type of electrostatic description of the chiral interaction 
can be combined with a Maier-Saupe mean-field treat- 
ment [s^ - lisj , or with a bare hard-core model and treated 
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within the seminal theory of Onsager [43 - l46| . 

A drawback of the Goossens potential is that its simple 
scaling form precludes an explicit connection with the un- 
derlying microscopic (electrostatic) interactions involved. 
As a consequence, the potential is unsuitable for studying 
the sensitivity of the cholesteric pitch of charged rod-like 
colloids where the strength of the chiral interactions de- 
pends strongly on the configuration of the surface charges 
and the electrostatic screening (viz. salt concentration). 
In this paper, we propose a simple but explicit model for 
long-ranged chiral interactions based on an impenetrable 
cylindrical rod decorated with a helical coil of charges lo- 
cated on the rod surface. The local interactions between 
the charged helical segments is represented by a simple 
Yukawa potential in line with the Debye-Hiickel approxi- 
mation valid for weakly charged poly-electrolytes. A sub- 
sequent analysis of the chiral potential between a pair of 
helical Yukawa segment rods leads to a simple general ex- 
pression which relates the intrinsic twisting potential to 
the electrostatic screening, charge density and coil pitch 
of the helices. The overall form of the chiral Yukawa 
potential resembles the one proposed by Goossens [36j . 
albeit with a much weaker decay with respect to the in- 
terrod distance. The most important difference, however, 
is that the amplitude of the chiral Yukawa potential can 
be linked to the coil configuration, and the local electro- 
static potential. By mapping the helical charge distribu- 
tion onto an effective helical coil the present model could 
be interpreted as a simple prototype for a charged twisted 
rod, a model commonl y in voked to explain cholesteric or- 
ganization of cellulose j28l . l47j and chitin microfibers [l^ . 
Likewise, the a-helical structure of polypeptide molecules 
[2?! could also be conceived as a coiled rod on a coarse- 
grained level. 

In the second part of the analysis, the potential will be 
employed in a simple Onsager second virial theory to pre- 
dict the behavior of the macroscopic cholesteric pitch of 
Yukawa helices. The sensitivity of the pitch with respect 
to the coil configuration and amplitude of the electro- 
static interactions will be scrutinized in detail. Also con- 
sidered are the implications of mixing two helices species 
with different lengths and/or coil pitches on the isotropic- 
cholesteric phase behavior of mixed systems. 

The theoretical results are consistent with experimen- 
tal findings and correctly capture the response of the 
cholesteric pitch upon variation of the particle concentra- 
tion and salt content. The theory also unveils a subtle re- 
lationship between the internal conformation (coil pitch) 
of the helical rod and the handedness of the cholesteric 
phase. The sensitive relationship between the symmetry 
of the cholesteric phase and the charge pattern on the 
rod surface is consistent with experimental observations 
in filamentous virus systems such as fd [i^ and MIS [i^ 
as well as numerical calculations based on a more explicit 
poly-electrolyte site model [i^ . 

This paper is structured as follows. In Section II we 
introduce the chiral Yukawa model and derive a general 
expression for the chiral potential imparted by the heli- 
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FIG. 1: Cylinders decorated with a continuous helical charge 
distribution (indicated by the thread) with different coil 
pitches, k — TT (left), k = 2n (middle) and k = An (right). 



cal charge distribution. The implications of the potential 
for the macrosopic cholesteric pitch will be analyzed in 
detail in Section III where we shall focus first on monodis- 
perse systems of identical helices followed by a treatment 
of binary mixtures of Yukawa helices differing in length 
and/or coil pitch. Finally, Section IV will be devoted to 
a discussion followed by some concluding remarks. 



II. CHIRAL INTERACTION BETWEEN 
YUKAWA HELICES 

In this study we aim to quantify the chiral or twist- 
ing potential of a pair of charged helical colloidal rods 
starting from a continuum electrostatic model based on 
a screened Coulombic (or Yukawa) potential. Within lin- 
earized Poisson-Boltzmann theory the electrostatic inter- 
action between two charged point particles with equal 
charge zLZe in a dielectric solvent with relative permit- 
tivity Er is given by [lOl : 



with r the distance between the particles, Xb = 
/4:7TeoerkBT the Bjcrrum length {Xb = 0.7 nm for 
water at T = 298-ftr), £0 the vacuum permittivity, and 
K the Debye screening constant k = AtiXb{Z p + 2cs) 
which measures the extent of the electric double layer 
surrounding each particle. Here, p = N/V is the macro- 
ion density, Nav is Avogadro's number, and Cg refers to 
the concentration of added salt. 

Eq. ([T]) can be generalized in a straightforward manner 
to describe the interaction between two helical macro- 
ions by assuming that the total pair potential can be 
written as a summation over Yukawa sites residing on 
the helices (see Fig. [ij . In the limit of an infinite num- 
ber of sites per unit length the electrostatic potential of 
helical rods with variable lengths Li and L2, and com- 
mon diameter D at positions {ri ,r2} and oricntational 
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unit vectors {ui U2} is then given by a double contour 
integral over the longitudinal (pitch) axes of the helices: 

^[/Y(Ar;ui,U2;Vi,V'2) = ^i^2Ab J dh J dt2 

^exp[-.|r^-rf|] 
\vl-vl\ 

with Ar = r2 — ri the distance vector between the ccntrc- 
of-masses of the helices, and Zj is the total number of 
charges on species i. 

In order to describe the segment position rf along the 
contour of helix i it is expedient to introduce a particle- 
based Cartesian frame spanned by tree orthogonal unit 
vectors {ui,v,Wi} {i = 1,2), where v = ui x U2/|ui x 
U2I and Wi = Ui X V. Within the particle the segment 
positions can be parametrized as follows: 

ri = ri + -2-iiUi + 2" 

X {cos(fciii + ^/'i)v + sin(A:iii + ?/'i)wi} 
L2 D 

1"2 = Y2 + —t2U2 + — 

X {cos(fc2<2 + ^2)v + sin(A:2i2 + i'2)^2} (3) 



in terms of the dimensionless contour variables ti (—1 < 
ti < 1) and coil pitch ki = 27rLi/p|^|. with the 
distance along the longitudinal rod axis over which the 
Yukawa helix makes a full revolution. With this con- 
vention fci > corresponds to a right-handed helix and 
fci < to a left-handed one. 

In this study we will focus on the general case of a bi- 
nary mixture of helical species with a different coil pitch 
(in terms of sign and/or magnitude), i.e. fci 7^ ^2. Since 
a helix is not invariant with respect to rotations about 
the longitudinal pitch axis we need to define an addi- 
tional orientational unit vector _L Uj to account for its 
azimuthal degress of freedom. Consequently, the interac- 
tion potential Eq. ([5]) must depend on a set of internal 
azimuthal angles ipi {0 < i^i < 27r), defined such that 
cos ipi = Wi • Gj . Using Eq. ([3]) and employing orthogonal- 
ity of the unit vectors allows us to derive an expression for 
the norm of the segment-segment distance. Let us nor- 
malize the centre-of-mass distance in units of rod length 
of species 1, Li, so that Ar — > Ar/Li. If we further 
assume the helices to be sufficiently slender {D/Li -C 1) 
we may Taylor expand the potential up to leading order 
in D/U 



/3f7y(Ar;ui,U2; V'l, V'2) 



/3f/^°'(Ar;ui,U2)-ZiZ2 



Xb 
D 



1) I "^^^I '^^^^y^^^'> 



X {[cos(fc2t2 + -02) - cos(fciii + 7/11)] (v • Ar) 

+ [sin(A:2i2 + ^2)(w2 • Ar) - sin(A:iii 4- V'i)(wi • Ar)] 

-(1/2) [ti sin(/c2i2 + ^/'2)(W2 • ui) + H2 sin(fci<i + Vi)(wi ■ U2)]} + 0[{D/Lf] (4) 

I 



in terms of the length ratio I ~ L2/ Li and Uy the deriva- 
tive of the Yukawa potential with respect to distance: 



[/{.(Af) 



exp[— K^iAf] / kLi 1 
2A^ \^ ^ Af2 



(5) 



and A7~ the distance between the positions ti and t2 along 
the pitch axis of the helices (in units Li): 



Af 



,~2 



Ar^ + [£<2(u2 • Ar) - ii(ui ■ Ar)] 
+ ][tl + ftj - 2£tit2{ui ■ U2)] 



(6) 



The reference potential fSUy in Eq. (j4]) represents the pair 
potential between two screened line charges of length Li 
and L2 interacting via the Yukawa potential (5l| : 

/3C/f (Ar;ui,U2) = ^1^2 (^^) f^dt, j'^dt2 

exp[— KLiAf] 



Ar 



Since this potential is strictly achiral we will disregard 
it in the sequel of this paper. Focusing now on the sec- 
ond, chiral contribution in Eq. ^ and applying standard 
trigonometric manipulations we may recast it into the 
following form: 



/3C/^'HAr;ui,U2;^i,V2) = 



D 



Li 



dt2C/y(Af) 



X {^cos?/;! + Bcos^(J2 + Csin?/;i + X'sin?/'2} (8) 



(7) Henceforth, we will denote this potential by superscript 
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"c" . The coefficients are given by: 

A = — cos(fciti)(v • Ar) — sin(A:iti)(wi ■ Ar) 

£t2 

— — sm{kiti){<vi ■ U2) 
B = cos(fc2t2)(v-Ar)+sin(fc2t2)(w2- Ar) 

sin(fc2^2)(w2 • Ul) 
C — sin(fciii)(v • Ar) — cos(fciti)(wi • Ar) 

cos(fciti)(wi ■ U2) 
V ^ - sin(fc2<2)(v • Ar) + cos(fc2t2)(w2 • Ar) 

COs(fc2t2)(w2 • Ul) (9) 

The next step is to construct an angle-averaged potential 
PUy by carrying out a proper preaveraging over the 
internal azimuthal angles. By requiring the Helmholtz 
free energy of the angle-averaged potential to be equal to 
that of the full angle- dependent potential one can show 
that 153: 




(10) 

where the brackets denote a double integral over the 
internal angles (.)^ = {27:)~^ J^^ dip. A similar ex- 
pression can be obtained starting from a self-consistent 
Boltzmann-wcighted average of the chiral potential, i.e 
Lf = (J7exp[— /3J7])^^ 2/(6xp[— /3C/])^^ 2 and Taylor ex- 
panding for small PU. This will give the same result ex- 
cept for a trivial prefactor in the fluctuation term [s^ [s^] ■ 

The leading order contribution {PUy'^)^^^^ vanishes [111 
upon integrating over ipi so that we need to consider 
the next-order fluctuation term in Eq. (jlOp . Using the 
isotropic averages (cos^ -0)^ = (sin^ ip)^ — 1/2 the angle- 
averaged chiral potential becomes: 

/34^)(Ar;ui,U2) - '{^1^1 [^)\^)' 

x{a^ + + & ^V^^ (11) 

where 

A = j dti j dt2U{.{AP)A (12) 

and identical relations for the other coefficients. A close 
inspection of the coefficients Eq. ([9]) reveals that only 
those terms proportional to the pseudo-scalar v • Ar con- 
tribute to the chiral potential. These are given by prod- 
ucts involving the first and third terms in Eq. All 
other contributions are invariant under a parity trans- 
formation Ar — Ar which renders them irrelevant 



for the present analysis. If we use the standard rep- 
resentation for the triple product Ui x (ui x U2) = 
ui(ui •U2) — U2(ui - Ul) so that — wi -112 = W2 -ui = sin 7 
(with 7 the angle between the main axes of the helices) 
the chiral potential simplifies to: 

^e<'',A,;a.a,,.-i^f^|(i^)'(^)^ 

xJ-i2(Ar;fci,fc2)x (13) 

in terms of a pseudo-scalar x which changes sign under 
a parity transformation [s^] Ar — Ar: 

X = (ui X U2 • Ar) (14) 

with Af = Ar/|Ar| the centre-of-mass unit vector. The 
function J-12 depends on the distance and helix orienta- 
tions and is invariant under parity transformation: 

J'i2(Ar;/ci,fc2) = -(cos(fciii))ti2(^t2sin(fci<i))t,2 

-(cos(A:2t2))ti,2 (^1 sin(fc2t2))ti,2 
-h(sin(fciii))ti 2 {£t2 cos(fciti))ti 2 

-h(sin(fc2<2))ti,2 {h cos(fc2i2))tiil5) 

where the brackets are short-hand notation for the double 
contour integration over pUy: 

, = Ari/2 dh J'^ dt2C/y(Af) (16) 

Since the prefactor J-12 depends rather intricately on the 
centre-of-mass distance and orientations it is desirable to 
seek a simplified form. This can be achieved by ignoring 
the interactions involving the ends of the helix. Indeed, if 
we take the limit of Li/D — > 00 the second contribution 
in Eq. which embodies the interaction between the 
end of one rod with the main section of the other, be- 
comes vanishingly small. An equivalent approach would 
be to fix the centre-of-mass distance vector along the unit 
vector Ar = Arv. In either case, the segment- segment 
distance Eq. ([6]) simplifies to: 

Af2 ^ + l{tl+ ftj - 2£tit2 COS 7) (17) 

and the only angular dependence is contained in the 
angle 7 between the main axes of the helices. The 7- 
dependence of the chiral potential at fixed centre-of-mass 
distance has been plotted in Fig. [2] and reveals a strongly 
non-monotonic relation. In a concentrated cholesteric ne- 
matic phase rods are usually strongly aligned along the 
local director so that 7 will on average be very small. In 
the asymptotic limit of strong local orientational order 
the following scaling expression for the chiral potential 
can be justified: 

/34'^ 0C7-Fi2(Ar;fci,fc2)(v- Af) (18) 

where A?^ = [Ar^ + \{ti- ^ts)^]^^^- Looking at the local 
minima appearing in Fig. [5] it is obvious that the asymp- 
totic approximation has to be taken with some caution in 
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the dilute regime where the average 7 is no longer very 
small. A striking anomaly occurs for fc = 4 where the 
chiral torque {dUy^ / 8^)^=0 has an opposite sign com- 
pared to the other cases shown. Moreover, the local and 
global minimum occurring for fc = 4 correspond to oppo- 
site twist directions. This hints to a subtle relationship 
between the magnitude of the coil pitch and the sense 
of the cholesteric director field which will be explored in 
detail further on in this study. In the asymptotic approx- 
imation, the distance dependence of the chiral potential 
is embodied by the chiral amplitude 7^i2 which will be an- 
alyzed in the subsequent paragraphs where we will focus 
on identical and enantiomeric helices, respectively. 

A. Identical helices 
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FIG. 2: Reduced chiral potential J^(Ar)|x| between two heli- 
cal Yukawa rods as a function of the interrod angle 7 for fixed 
centre-of-mass distance Ar = O.IL and nL = 20. 



If the helices are identical, then Li = L2 = L {£ = 1), 
ki = k2 = k and Z\ = Z2 = Z . The double contour 
integration Eq. (|16|) is invariant under interchanging t\ f-> 
^2 and we can exploit this symmetry to simplify Eq. (|15p 
as follows: 

^(Ar;fc) = -2(cos(Hi))t,_,(i2sin(fcti))t,, 

+2(sin(fc<i)}t, , (t2 cos(Hi))t, , (19) 

The second contribution involves odd terms in t\ and ^2 
which vanish upon performing the double contour inte- 
gration [Eq. ([T6)) ]. The result is: 



T{I^r;k) = -2Ar^y" dh j dt2?7y (Af ) cos(fcti) 

^ (/ ^*^/ '^^2J/y(Af)t2sin(Hi)^(20) 

which still implicitly depends on the coil pitch k and De- 
bye screening length kL via Eq. ([5]) . The chiral potential 
between two identical helical Yukawa rods thus reads: 



/34'^(Ar;ui,U2) 
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As 
D 



^] H^r-k)x 



(21) 

This potential, shown in Fig.[3l is found to decay steeply 
with increasing rod centre-of-mass distance. Owing to 
the intractable double contour integrations, the distance- 
dependence of the potential is essentially non-algebraic 
and does not obey a simple power-law scaling with re- 
spect to Ar. At short distances the amplitudes appear 
to be mainly dependent upon the coil pitch k rather than 
the screening constant. 

It is worthwhile to compare the chiral potential 
Eq. (PT|) to the one proposed by Goossens [H, I3 which 
has been used frequently in literature. The Goossens po- 
tential emerges as the leading order contribution from a 
multipolar expansion of the interaction potential between 
two molecules each composed of an array of dipoles. This 
potential takes the following generic form: 

-U^[-^)\n,-n2)x (22) 
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FIG. 3: Amplitude of the chiral potential [Eq. (f20)) ] between 
two helical Yukawa rods as a function of centre-of-mass dis- 
tance Ar = Ar ■ v. 



in terms of a reference size a and amplitude C/q which 
determines the handedness of the chiral interaction. As 
can be gleaned from Fig. [3] the distance-dependence of 
the Goossens potential is much steeper compared to the 
decay found from Eq. (|20|) . We further remark that both 
chiral potentials are invariant with respect to an inversion 
of the rod orientation (ui — >■ — Ui, or equivalently 7 -^• 
7 ± tt). 



C/Gs(Ar; Ui,U2) 



B. Enantiomers: helices with opposite handedness 



Let us now consider two helices of identical length Li = 
L2 ~ L and charge Z\ = Z2 = Z but opposite pitch 
k\ = — fc2. The corresponding cross interaction between 
a right-handed and left-handed helix then follows from 
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Eq. (HID: 



J"i2(Ar; fci, -fci) 



-(cos(fciti))ti 2(^2 sin(/citi))tj_2 
-(cos(-fcit2))fi_2(ti sin(-fcit2))ti_2 
+ (sin(A:iti))tj 2 {t2 cos{kiti))t^ .^ 
+ (sin(-fcit2))ti,2 (*i cos(-fcit2))ti 2 
(23) 



irrespective of the magnitude of the coil pitch. The impli- 
cation of this result is that an equimolar binary mixture 
of helices with equal shape but opposite handedness (a 
so-called 'racemic' mixture) does not exhibit cholcstcric 
nematic order. If the rods have different lengths {£ ^ 1) 
they no longer form an enantiomeric pair and the chiral 
interaction will generally be non-zero. 



III. PREDICTION OF THE CHOLESTERIC 
PITCH 



In order to study the implications of the internal heli- 
cal structure of the rods on the (macroscopic) cholcstcric 
order we need a statistical theory which is able to make 
a connection between the pair interaction of the parti- 
cles and the local equilibrium orientational distribution 
/(u • fi) and cholcstcric pitch. Such a theory can be de- 
vised starting from the classical Onsager theory [i^] for 
infinitely thin hard rods which exhibit common uniaxial 
nematic order. The theory has been generalized by Stra- 
ley [55I [scj for aligned fluids with weakly non- uniform di- 
rector fields, such as a cholesteric liquid crystal, in which 
case elastic contributions must be incorporated into the 
free energy. The Helmholtz free energy density F/V of 
a binary mixture of slender chiral rods in a cholcstcric 
phase of volume V takes the following form: 



V 



= p{ln pV - 1) + pJ2 J dufi{u)\n[x^4:7^f^{u)] 



(24) 



In terms of the thermal energy {3^^ ~ , overall parti- 
cle density p = N/V. V = Y\ - V^^ is a weighted product 
of the thermal volume Vi of particle i which includes con- 
tributions arising from the rotational momenta. The first 
two terms in the free energy denote the ideal translational 
and orientational entropy of the system while the last one 
represents the excess free energy which accounts for the 
interactions between the rods on the approximate second- 
virial level. The latter consists of three contributions. 
The first, Kq, involves a spatial and orientational aver- 
age of the Mayer function /m = exp[— /3J7] — 1 weighted 
by the orientational distribution functions (ODF) fi{u) 



of the respective species 



C^Ui/,(Ui) J dU2fj{u2) 

J dAr/7,(Ar;ui,U2) 



(25) 



For hard anisomctric bodies, the spatial integration leads 
to the excluded volume w^^^[(ui,U2) between particles of 
species i and j. The second and third contributions in 
the excess free energy represent the change of free en- 
ergy due to the twist deformation of the director field. 
The strength of this deformation is measured by the 
cholesteric pitch q = 27r/p, with p ^ L ^ pint the pitch 
distance associated with the helical director field. Since 
the theory is only valid for long- wavelength distortions of 
the director field it is required that g ^ 1. The torque- 
field contribution, proportional to q arises from the 
chiral torque imparted by the chirality of the particles 
and leads to a reduction of the free energy. Opposing 
this, there is an elastic response counteracting the de- 
formation of the director field. The corresponding free 
energy penalty K2q^ is proportional to the twist elastic 
constant K2 ■ Similar to Kq, the torque-field and twist 
elastic constants are given by spatio-angular averages of 
the Mayer function, albeit in a more complicated way 

= -V / d-i^ihi-ixi) I du2/,(u2)(u2-y) 



j dAr(Az)/li(Ar;ui,U2) 



dui/,(ui)(ui -y) / rfu2/j(u2)(u2 -y) 



(26) 



j dAr{Azff%iAr;u,,n2) 



where / represents the derivative of the ODF with re- 
spect to its argument. In arriving at Eq. ([26)) we have 
implicitly fixed the pitch direction along the z-direction 
of the laboratory frame with the local nematic director 
n(z) pointing parallel to the x-axis. 

Let us now consider a model binary mixture of hard 
rods with different lengths (£ 7^ 1) decorated with a weak 
chiral potential of the form proposed in the previous sec- 
tion [Eq. p^ ]. The twist elastic constant is not affected 
by the weak chiral potential, but only by the achiral hard 
cylindrical backbone and electrostatic reference poten- 
tial. If we neglect the latter contributions for the time 
being [gHI the spatial integrals appearing in Eq. (f25|) and 
Eq. (j26p reduce to integrals over the excluded volume 
manifold of two thin hard cylinders weighted over pow- 
ers of the pitch distance variable Az: 



M^(ui,U2)= / Ar(Az)", 



n = 0,2 



(27) 



These quantities have been calculated in general form for 
hard spherocylinders in Ref. |4^ Here, we only need the 
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leading order contributions for large Li/D: 



coil of species i: 



M^J(ui,U2) = t;;4^,(ui,U2) =2L,LjD|sin7| 



r2 
4 



+-^(u2-zf + i?^(v.zf (28) 



Due to symmetry reasons the torque-field constant Ki 
depends only the chiral part of the potential and not 
on the achiral reference part. For weak chiral potentials 
considered here it is justified to approximate /m k, — jSU, 
analogous to a van der Waals perturbation approxima- 
tion generalized to liquid crystals [i^ll^,!!!] ■ The spatial 
integration pertaining to Mi then becomes: 



Afi^^(ui,U2) 



Ar(Az);3i7^''Vr;ui,U2) (29) 



where the spatial integral runs over the space complemen- 
tary to the excluded volume of the particles. By exploit- 
ing the cylindrical symmetry of the cholesteric system one 
can parametrize the distance vector in terms of cylindri- 
cal coordinates so that Ar = r sin i^x -|- r cos Qy + Azz. 
With this, one can write [67| : 

M^^(ui,U2) = iz2z2A|i?2y^(fc„fc,)(ui X U2 ■ z) (30) 

where W represents a spatial integral over the amplitude 
of the chiral potential for a given pitch ki of the Yukawa 



W{k,,kj) = -An drr dAz{Az)^j:,j{Ar;h,kj) 

Jo JD/Li 

(31) 

with Ar = (r^ + Az^)^/^ the centre-of-mass distance 
parametrized in terms of cylindrical coordinates Ar and 
r (both in units of Li). 

The next step is to perform a double orientational 
averages of the moment contributions according to 
Eq. ((26|) . It is expedient to adopt the Gaussian approx- 
imation, in which the ODF is represented by /i(ui) oc 
exp[— ai(ui • n)^/2] in terms of a single variational pa- 
rameter ai whose equilibrium value is required to mini- 
mize the total free energy. If the local nematic order is 
strong enough (a^ 3> 1) the orientational averages can 
be estimated analytically by means of an asymptotic ex- 
pansion for small inter-rod angles. This procedure has 
been outlined in detail in Refs. ll^. The result for the 
nematic reference contribution i^Q-' reads (up to leading 
order in a;): 



!3Ki' ^ p'L,L,D{{j)) 



p'ULjD 



1/2 

^^ ^ "J \ (32) 



where the double brackets denote Gaussian orientational 
averages, specified in the Appendix. Similarly, one can 
derive for the twist elastic constant: 



p 

192 



UL,Da.,a, {l\ {{{^9^0^ + 0]))) - ((7^0?))) + {{{j0^{0^ + 0^))) - ((7^0|)))] 

4Q!iQ!j)\ 



(33) 



192 



1/2 1/2, 
OL- a.' (a. 



a. 



)3/2 



Here we have ignored the last term in M2 [Eq. ([28)) ] 
which is negligible for slender rods. Finally, the torque- 
field contribution reads: 



(34) 



With this, the free energy is fully specified. Minimization 
with respect to q yields for the equilibrium pitch: 



q 



(35) 



For weak pitches [q <C 1), the local nematic order is 
only slightly affected by the twist director field and the 
equilibrium values for ai depend only on the free energy 



of the nematic reference phase {q — 0). Minimization 
with respect to ai and rearranging terms leads to the 
following set of coupled equations: 



ai 



Q 



Q 



-lN-1/2 



2i/V£(l + Q)-V^ + :r2^^ 
a;i-h2i/2a;2£(i + Q-i)-i/2 



(36) 



in terms of the ratio Q ~ a2/ai, dimcnsionlcss concen- 
tration c = NLiD/V and length ratio £ = L2/L1. These 
equations cannot be solved analytically but the solutions 
are easily obtained by iteration. It is important to note 
that both Qfi and a2 increase quadratically with the con- 
centration c since their ratio Q only depends on the mole 
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fractions x, 



A. Monodisperse systems 

For pure systems of infinitely thin helices, the twist 
elastic constant behaves asymptotically as: 



10 



PK2D 



14c 
192 



(37) 



as found by Odijk Ut 
tain: 

qL = 



For the cholesteric pitch we ob- 



PKiDL 12 



liK2D 



7 



c£c 



(38) 



where £'c(Ar) is a dimensionless parameter pertaining to 
the chiral potential between rods. It combines the micro- 
scopic characteristics of the helical rod such as the aspect 
ratio, surface charge and coil pitch [viz. Eq. ([50))]: 



(39) 



Since the electrostatic screening of the charge-mediated 
chiral interactions depends on the rod concentration the 
variation of the pitch with concentration is strictly non- 
linear, as shown in Fig. |4l To simplify matters we may 
state that typically Xb/D ~ 0(10""'^) for rod-like col- 
loids in water. The rod charge Z is expected to be lin- 
early proportional to the rod length L. Let us further 
introduce a charge density Cc, defined as the number 
of unit charges per unit length, so that Z ~ UcL and 
Z'^{D/LY ~ 0{L/D) (assuming the charge density UcD 
to be of order unity). We remark that the qL values 
shown in Fig. 2] correspond to pitch distances of 10-100 
rod lengths, a range which is commonly found in experi- 
ments |16l.[l7|. 

Fig. 0] shows that the pitch is a monotonically increas- 
ing function of the rod concentration with the magnitude 
depending strongly on the salt concentration. At high Cg 
the electrostatic repulsion between the Yukawa sites on 
the rods is strongly screened and the resulting chiral in- 
teraction will be attenuated significantly. The screening 
effect is common in dispersions of fd and DNA and sup- 
ports the idea that chiral interactions are primarily trans- 
mitted via long-range electrostatic interactions [3, . 
We remark that for short-fragment (146 base pair) DNA 
the opposite trend is observed [l^. There, the cholesteric 
pitch is found to increase with respect to salt concentra- 
tion, i.e. charge screening leads to stronger chiral interac- 
tions between the DNA chains. This trend is most likely 
explained by the steric effect associated with the double 
helical backbone which becomes more pronounced as the 
charged phosphate groups residing on the backbone be- 
come increasingly screened. Clearly, the excluded volume 
of the helical grooves (neglected in this study) must be 
included explicitly to give a proper account of both steric 
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FIG. 4; Variation of the cholesteric pitch qL with concen- 
tration c for right-handed Yukawa helices with fc = 1.4 and 
L/D = 50. The charge density is (TcD = 1 and the reduced 
Bjerrum length Xb/D = 0.1. 



and electrostatic contributions to microscopic chirality in 
DNA 



35|. 



We remark that the effect of the temperature on the 
chiral strength is rather trivial in our model. At high 
temperatures the chiral dispersion potential will become 
increasingly less important than the achiral hard-core po- 
tential associated with the cylindricasl backbone and a 
reduction of the pitch is expected. This effect has been 
observed in fd virus rods [13, Ell . 

The symmetry of the chiral interaction is entirely gov- 
erned by the spatial integral over the chiral potential 
which depends rather intricately on the coil pitch of the 
Yukawa helix and the Debye screening length. Fig. [5] 
shows the typical behavior of the intrinsic chiral strength 
for two helical rods with i/D = 50 as a function of the 
coil pitch. For small coil pitches, i.e. weakly coiled heli- 
cal rods, the handedness of helical director field is com- 
mensurate with that of the Yukawa coil, e.g right-handed 
helices form a right-handed cholesteric phase. The max- 
imum twisting effect is obtained for k = 1.35. This value 
corresponds to a coil pitch distance of pint ^ 5L which 
implies that only a marginal degree of coiling is required 
to minimize the cholesteric pitch. 

A further increase of k leads to a reduction of the chiral 
strength and a sign change at fc = 3.26 corresponding to 
a sense inversion of cholesteric helix. For |fc| > 3.26 the 
handedness of the cholesteric helix is no longer commen- 
surate with that of the Yukawa helix so a left-handed 
cholesteric state is obtained from right-handed helices 
and vice versa. Increasing fc even further reveals a oscil- 
latory relation between the microscopic (coil) and macro- 
scopic (cholesteric) pitches. The inversion points depend 
primarily on the coil pitch and are not affected by the 
electrostatic screening kD. At high fc, the twisting po- 
tential strongly decays and vanishes asymptotically for 
tightly coiled Yukawa rods (fc — )■ 00) where the effective 
'width' of the helical grooves becomes negligibly small. 
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FIG. 5: Spatially integrated chiral potential >V(fc) [Eq. jSTJ] 
versus coil pitch k — 2-KL/pint for identical Yukawa helices 
with a right-handed (fc > 0) or left-handed {k < 0) symmetry 
for L/D = 50. For W > a right-handed cholesteric phase 
is formed, whereas W < corresponds to a left-handed one. 
Zero-points indicate a sense inversion of the cholesteric helix 
at A; = 3.25,4.67, 6.51,8.00, ... 



B. Binary mixtures of helical rods with equal 
lengths = 1) 

For a mixture of rods with equal length but different 
coil pitches (fci ^ ^2) and/or surfaces charges {Zi ^ 
Z2), the situation is comparable to the monodisperse case 
since the twist elastic constant is independent of these 
properties. Furthermore, we have ai ~ a2 = o: and the 
cholesteric pitch is given by a form analogous to Eq. 



qL 



12c 

T 



(40) 



in terms of an effective chiral strength given by a simple 
mole fraction average of the different pair contributions: 



X'lXj 



(41) 



with the generalized version of Eq. ([5^ : 

2 / \ 3 



£'J 



Xb 
D 



ZfZ] 



W{k,,k,) (42) 



C. Binary mixtures of helical rods with unequal 
lengths / 1) 

In case of a binary mixture of rods with different 
lengths things are more complicated and most thermo- 
dynamic properties such as the cholesteric pitch can only 
be assessed numerically. Let us first investigate the ef- 
fect of a weak chiral potential on the phase behavior of 
rod mixtures with different length ratios. To that end 
we must first consider the free energy of a binary mix- 
ture of infinitely thin hard rods with length ratio i in the 



ncmatic phase. Within the Gaussian approximation it is 
given by [s^ [gO] : 



N 



In c -f Xi [In XiCii — 1] 

i 



+^EE 



Li 



1/2 



^ In c + Xi [In 

i 

X \xl + 2^^'^xxX2l{\ + 0"^)^/^ + x\l^Q-\il 
In the isotropic phase the free energy simplifies to: 

^ In c + Xi In Xi + c— [xf + 2xiX2( + x'^i'^] 

(44) 

In the cholesteric phase we must take into account the 
change of free energy associated with the weak helical 
distortion of director-field as discussed in the preceding 
Section: 



chol 



N 



N 



7EE 



XiXj 



(45) 

which yields upon minimization with respect to the 
cholesteric pitch q = 2iTLi/p: 



N 



N 



1 {LiJ2jX^X3Kl^ 



(46) 



Wc reiterate that this expression is only applicable in the 
regime 5 <C 1 where the helical distortion of the director 
field does not interfere with the local orientational order. 
Since the twist elastic constants are positive, the second 
term in the free energy must be negative which implies 
that a small degree of chirality always leads to a reduction 
of the free energy of the system. In explicit form, the 
twist elastic constants for the pure components are given 
by: 



K. 
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K. 



22 



(3K'2'D 



If. 

192 Vai 
192 \ai 



1/2 



1/2 



-1/2 



K, 



12 



192 V ai 



1/2 



g{Q) 



(47) 



with 



.(.^ gi/2 (3Q^ + 4Q)^+(3 + 4g)^=^ f48^ 
^(^^"^ 7QV2(i + Q)3/2 (48) 

Similarly, one can produce for the torque-field contribu- 
tions: 
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With this, the cholesteric free energy of a binary mixture 
with 7^ 1 can be rewritten as: 



~/V TV "28 



with: 



G{x,,Q) 



xi+2V^X2^(l + Q-i)-i/^ 

x2+2xiX25((3)+Xp4Q-l/2 



(49) 



(50) 



The corresponding cholesteric pitch takes the following 
form: 



, 2^Li 12c- 

q = —£cG{xi,Q) 

V 7 



(51) 



As required, this expression reduces to the forms pro- 
posed in the previous subsections [Eq. (j3Hl) and Eq. (HOI) ] 
if we substitute £ = 1 and Q = 1. Recalling that 
Q = Q{x) wc conclude that the concentration depen- 
dence of the cholesteric pitch of the mixture is identical 
to that of the monodisperse system. The behavior of the 
pitch with respect to the mole fraction 12 = 1 — xi on 
the other hand is expected to be quite rich and this will 
be scrutinized in detail further on. 

Let us now consider a binary mixtures of helical rods 
with equal charge density <Tc so that the mole-fraction 
averaged chiral parameter Be can be written in compact 
form as: 



£c = eo {x\eii + 2.Tia;2£^ei2 



£22 



(52) 



in terms of an amplitude Sq and phase factor Si 



£0 



d) \d 



(53) 



which needs to be specified for the mixture of interest. 

Let us first consider the case £ = 3 with ki = —k2 = 
1.4, i.e. a mixture of short right-handed helices mixed 
with long left-handed ones with equal pitch magnitude. 
To simplify matters we shall consider dispersions with 
excess added salt such that the screening constant is 
fixed to kLi = 50, independent of the rod concentra- 
tion. The phase diagrams of the chiral systems and the 
corresponding hard rod reference system are shown in 
Fig.ini Althought the isotropic-nematic transition is only 
marginally affected by chirality, the cholesteric phase 
show a dcmixing into two coexisting cholesteric fractions 
with opposite handedness. The dcmixing region is closed 
off by a lower critical point located at = indicating 
that the cholesteric pitch vanishes at the critical point. 
Note that the achiral mixture do not show this dcmix- 
ing at this particular length ratio. Within the Gaussian 
approximation (60| such a demixing only occurs above a 
critical length ratio i > 3.167. Moreover, the nematic- 
ncmatic binodals do not meet in a lower critical point 
but merge with the isotropic-nematic ones to produce 



an isotropic-nematic-nematic triphasic equilibrium, irre- 
spective of the length ratio i. The cholesteric-cholesteric 
demixing is driven primarily by the small chiral disper- 
sion contribution to the rod interaction potential and 
does not arise from an interplay of the various entropic 
contributions (associated with mixing, free-volume and 
orientational order) such as for hard rod mixtures. A sim- 
ilar demixing is observed for a mixtures of right-handed 
helices with fci fc2 = 1.4 at £ = 3 and amplitude 
eo = 0.1, albeit at a much higher osmotic pressure. It 
is worth noting that for this case the cholesteric pitch 
does not reduce to zero at the critical point but attains a 
finite positive value. For the case £ = I, i.e. helices with 
equal lengths, no demixing of the cholesteric state was 
found. This shows that a cholesteric mixture of left- and 
righthanded helices can only demix if the length ratio 
between the species is considerable. 

Fig.[7|shows the behavior of the cholesteric pitch across 
the range of mole fractions. For the first mixture, a quasi- 
linear reduction of the pitch is observed accompanied by 
a change of handedness upon increasing the mole frac- 
tion of the left-handed 'dopant' helices. The zero-point 
corresponding to a vanishing pitch is virtually (but not 
strictly) independent of the osmotic pressure of the sus- 
pension. For the second mixture the trend is completely 
different and features a non-monotonic increase of the 
pitch with mole fraction. A similar non-monotonic trend 
is observed for the pitch of a cholesteric phase in coexis- 
tence with an isotropic phase as a function of pressure. 



IV. DISCUSSION AND CONCLUSIONS 

It is instructive to compare the present Yukawa-type 
chiral potential with a much simpler model potential em- 
ployed in a previous study Ref . l46l There, we have used a 
square well (SW) chiral potential based on the Goossens 
form Eq. 



;7sw(Ar; Ui, U2) = -eswH{Ar - A)(ui • U2)x (54) 

with H{Ar — A) a Heaviside step fimction. The main pa- 
rameters characterizing the chiral interaction arc the SW 
range A > {1 + D/L) and depth esw- The main advan- 
tage of using the SW form is it that it renders the spa- 
tial integrations over the potential analytically tractable. 
The result for the torque field constant is very simple: 
PKiD^ ~ (7rcV6)/3£swA''. Comparing this with Eq. ^ 
and setting the reduced SW range A to unity allows us 
to make an explicit link between the SW depth and the 
microscopic parameters pertaining to the electrostatic in- 
teractions between the helices: 

The justification for this relation lies in the notion that 
most thermodynamic properties are governed by the spa- 
tially integrated pair potential rather than the bare one. 
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FIG. 6: (a) Isotropic-cholesteric nematic phase diagram for a mixture of Yukawa helices with length ratio £ = 3 consisting 
of helices with opposite handedness (fci = — — 1.4). Plotted is the osmotic pressure P = piiLlD versus the mole fraction 
of the left-handed species x = X2- (b) High pressure region featuring a demixing of the cholesteric phase. A critical point is 
indicated by the dot. (c) Cholesteric pitch of coexisting left- and right-handed cholesteric phases versus pressure, (d)-(f) Same 
results for a mixture of helices with equal handedness (fci — k2 = 1.4) and £ = 3. 




FIG. 7: Variation of the cholesteric pitch versus mole fraction 
at constant pressure, (a) Helices with opposite handedness 
(fci = — = 1.4, ^ = 3). (b) Helices with equal handed- 
ness (fci — k2 — 1.4, £ — 3). (c) Behavior of the pitch as a 
function of the isotropic-cholesteric coexistence pressure for 
ki = -k2 = 1.4 and £ = 3. (d) Same for fci = fca = 1.4, £ = 3. 



A prominent example is the classical van der Waals model 
for fluids whose universal nature stems from the fact 
that any arbitrary (but weakly) attractive potential can 
be mapped onto a single integrated van der Waals en- 
ergy which, along with the excluded-volume contribu- 
tion, fully determines the equation of state and hence 
the thermodynamics of the fluid state. 

In summary, wc have proposed a helical Yukawa seg- 



ment model as a course-grained model in an effort to 
quantify the chiral interaction between stiff chiral poly- 
electrolytes. Our chiral potential strongly resembles the 
classic Goossens potential [s^, which has been used al- 
most exclusively to describe long-ranged chiral disper- 
sion forces. Contrary to the Goossens form, our poten- 
tial provides an explicit reference to the microscopic and 
electrostatic properties of the rods. Combining the po- 
tential with a simple second-virial theory allows us to 
study the structure and thermodynamic stability of the 
cholesteric state as a function of rod density, degree of 
coiling, and electrostatic screening. While the magni- 
tude of the cholesteric pitch depends strongly on the 
rod concentration and concentration of added salt, the 
handedness of the phase is governed mainly by the pitch 
of the Yukawa coil. The symmetry of the cholesteric 
phase need not be equivalent to that of the individual 
helices but may be different, depending on the precise 
value of the coil pitch. Within a certain interval of the 
coil pitch, right-handed Yukawa coils may generate left- 
handed cholesteric order and vice versa. The antagonis- 
tic effect of charge-mediated chiral interactions is consis- 
tent with experimental observations in Mi 5 virus systems 
[igt and various types of DNA [s^] where left-handed 
cholesteric phases are formed from right-handed helical 
polyclectrolyte conformations. Small variations in the 
shape of the helical coil, induced by e.g. a change of tem- 
perature, may lead to a sense inversion of the cholesteric 
helix. Such an inversion has been found in thermotropic 
(solvent free) polypeptides [6l[ and cellulose derivatives 
[62 | , and in mixtures of right-handed cholesterol chloride 
and left-handed cholesterol myristate . 

Our model also predicts that a very small degree of 
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microscopic hclicity is required to maximize the twisting 
potential of the Yukawa hehces. The optimum is reached 
when the pitch distance of the coils equals about 5 times 
the rod length. Mixing stiff helical rods with sufficiently 
different lengths may lead to a demixing of the cliolesteric 
phase at high pressures. The demixing region closes off 
at a critical point upon lowering the osmotic pressure. 

The present model could be interpreted as a simple 
prototype for complex biomacromolecules such as DNA 
and fd which are characterized by a helical distribution of 
charged surface groups. Other lyotropic cholesteric sys- 
tems, such as cellulose and chitin microfibers in solution 
could also be conceived as charged rods with a twisted 
charge distribution (T^]. Small changes in the internal 
twist of the fibrils could be induced (e.g. by applying 
an external field or varying the temperature) in order to 
tunc the handedness of the cholesteric phase. This could 
be of importance for the use of chiral nanocrystals in op- 
tical switching devices and nanocomposites. We remark 
that a more accurate description of the pitch sensitivity, 
particularly for DNA systems, could be formulated by 
accounting for the steric contributions associated with 
the helical backbone of the chains as well as the influ- 
ence of chain flexibility. This could open up a route to- 
wards understanding the unusual behaviour of the pitch 
versus particle and salt concentration as encountered in 
DNA [6|,|l| 

using simple coarse-grained models. Further- 
more, a quantitative comparison of the predicted pitch 
distances with experimental systems should be possible 
for rigid, slender polyelectrolytcs with a well-defined sur- 
face charge and internal pitch. Most chiral systems stud- 
ied to date however do not fulfill these criteria which 
makes it difficult to put our predictions to a quantitative 
test. At present, our theory should therefore be consid- 
ered as a mere qualitative guideline. 
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Appendix: Gaussian averages 

The Gaussian averages required for the calculation of 
the twist elastic constant arc taken from Rcf. [6^ We 
quote them here: 
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